Dark Matter and Gauge Coupling Unification in A 
Supersymmetry Model with Vector-like Matter 

Chun Liu and Jia-Shu Lu 

State Key Laboratory of Theoretical Physics and 
Kavli Institute for Theoretical Physics China, 
Institute of Theoretical Physics, Chinese Academy of Sciences, 
P.O. Box 2735, Beijing 100190, ChinJA 
(Dated: May 2, 2013) 

Abstract 

WIMP dark matter and gauge coupling unification are considered in an R-parity violating MSSM 
with vector-like matter. Dark matter is contained in an additional vector-like SU(2)^ doublet which 
possesses a new U(l) gauge symmetry. The Higgs fields are extended to be in a 5 © 5 representation 
of SU(5). The stability of dark matter is a result of gauge symmetries, and the mass of the 
dark matter particle is between (1.1-1.5) TeV. Dark matter has a very small cross section with 
nucleis, thus the model is consistent with current dark matter direct detection experiments such 
as XenonlOO. The model also predicts new charged and colored particles to be observed at LHC. 
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I. INTRODUCTION 

Dark matter is one of the most important problems in elementary physics It plays a 
very important role in understanding some astro-physical observation. If dark matter par- 
ticles are thermally produced in the early Universe, an attractive scenario appears, namely 
;hat the dark matter particle has a hundreds GeV mass with typically weak interaction 
2)]. The weakly interacting massive particle (WIMP) scenario of dark matter is interest- 
ing for particle physics. In high energy physics experiments searching for the electro-weak 
symmetry breaking (EWSB) mechanism, WIMP should be found in the near future. 



lowing main stream 
j| . Supersymmetry 



The particle physics beyond the Standard Model (SM) has the fo 
logic. SM gauge interactions unify at a high scale (~ 10 16 ) GeV 
(SUSY) 4[ is then required to stabilize the Higgs mass. The minimal SUSY extension of 
SM (MSSM) which necessarily involves two Higgs doublets makes the idea of the grand 
unification theory (GU r 



the lightest neutralino 



') more meaningful due to LEP data [5]. WIMP dark matter is right 
2J by further assuming R-parity conservation, . 
We will work in the SUSY paradigm without assuming R-parity conservation. R-parity 
conservation is usually adopted to avoid rapid decays of the proton, however, it lacks of 
motivation from first principles. Instead, we can only assume baryon number conservation 
which is also phenomenologically viable. R-parity violation makes things more complicated. 
The lightest neutralino is no longer stable. To have dark matter, new particles are needed 
then. To keep dark matter WIMPs, the simplest realization is to introduce a vector-like 
SU(2)l doublet, like the two Higgs doublets in MSSM, with a weak scale mass. We take one 
fermionic neutral component of the new doublet superfields as the dark matter particle. Its 
stability requires a new symmetry which is a gauge symmetry. The necessity of a new gauge 
symmetry, instead of a discrete symmetry, lies in the following point. The elastic scattering 
of the dark matter particle via a Z boson exchange should be suppressed. In our case, the 
dark matter particle is not really a Dirac one, because spontaneous breaking of the new 
gauge symmetry splits the neutral particle spectrum. As a result, dark matter and the Z 
boson interaction is almost inelastic. 

Now we consider how to let the new dark matter particle content compatible with gauge 
coupling unification. Because one vector-like doublet is introduced, in addition to the par- 
ticles of MSSM, gauge coupling constants would no longer unify in this case. GUT relation 
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of gauge coupling constants could be restored if some colored particles are further added so 
that all the new particles form complete representations of GUT (], 7|. Alternatively, we can 
keep the dark matter sector as simple as possible, namely it is just SUSY two-Higgs-doublet 
alike, but we attach the new colored particles with the two Higgs doublets. In an effort of 
extending MSSM, we proposed that Higgses are understood as sleptons of an extra vector- 
like generation [8]. In that model, GUT was lost, and there was no dark matter candidates. 
It is interesting to note that whence we consider the above discussed dark matter scenario, 
GUT relation of the gauge coupling constants can be restored. This will make the whole 
low energy SUSY model more meaningful. 

However, careful consideration about running gauge coupling constants tells us that 
things are not that easy and straightforward. The particle content of the extra vector- 

n 

like generation [8J needs to be reduced. To avoid Landau poles of the coupling constants, 
new particles at the TeV scale cannot be that many. Instead of a TeV vector-like generation 
which would have both a 5 © 5 and a 10 © 10 representations, we are only allowed to have 
a TeV 5 © 5 in which the two Higgs doublets are contained. Therefore, as far as the Higgs 
content is concerned, we return back to the ordinary GUT. Here the new point is that the 5 
representation will mix with ordinary three generation fermions, and the triplet Higgs mass 
is around TeV. We impose by hand at the moment baryon number conservation. 

In the next section, we present the model which includes the dark matter sector. In Sect. 
Ill, dark matter properties and collider phenomenology are analyzed. Discussions and the 
summary are given in the last two sections. 



II. MODEL 



Let us first look at the SM relevant sector. Within the framework of SUSY and 
SU(3) c xSU(2)ixU(l)y gauge symmetry, the particle content is extended in a way that 
the two Higgs doublets are contained in a SU(5) 5 © 5 representation. These vector-like 
particles have masses of (100-1000) GeV. Because of R-parity violation there is no dark 
matter particle. 

To include the dark matter sector, we note that all the particles already introduced are 
in full representations of SU(5) GUT. For the GUT purpose, the WIMP sector should be 
exactly composed of a pair of SU(2)^ doublets with opposite U(l)y charges. This is the 
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reason to take them SUSY two-Higgs-doublets-like. To avoid extra degrees of freedom which 
may violate GUT relation, the new interaction among the WIMPs should be another Abelian 
gauge interaction U(l) n . The U(l) n charge can be arranged so that after U(l) n breaking, 
an unbroken Z 2 symmetry will remain. This Z 2 symmetry makes dark matter stable. 

The model is SUSY SU(2) i xU(l)yxSU(3) c xU(l) n gauge symmetric one with baryon 
number conservation. The particle content is given below with their quantum numbers 
under above gauge symmetries and global baryon numbers, 

L m (2, -1, 1, 0, 0) , E<(1, 2, 1, 0, 0) , Qi{2, §, 3, 0, §) , Uf (1, -|, 3, 0, -|) , D c m (l, §, 3, 0, -|) , 

# u (2,l,l, 0,0), L^(l,-§,3,0,±), 

(1) 

and 

X i(2, -1, 1, -1, 0) , X2 (2, 1, 1, 1, 0) , 0, 1, -2, 0) , 2 (1, 0, 1, 2, 0) , X(l, 0, 1, 0, 0, 0) . (2) 

Field notation is conventional, like L stands for lepton doublets, E c for lepton singlets, Q 
for quark doublets, U c - D c for quark singlets and H u for the up-type Higgs doublet. In Eq. 
(PQ), m = 1 — 4 and i = 1 — 3. One combination of L m 's makes the down- type Higgs. In 
addition to the ordinary three generations, there is a vector-like 5-representation (L4, D\ 
and H u , D C H ). The dark matter sector is given in Eq. ((2]). Two doublets X1.2 contain the 
WIMP, and 0i i2 are the U(l) n breaking Higgs. 



A. SM relevant part 

The superpotential of our sector containing the SM part can be written down. Instead 
of R-parity, baryon number conservation is assumed, 

W = fi m L m H u + ^D c m D c H + \ mni L m L n El + X' imn QiL m D^ + y tj Q t H u U< + y t ^D c H U^ , (3) 

where /i m 's are mass parameters, A^, y and y's coefficients. By redefining the down-type 
Higgs and the fourth down-quark field, 

Hd = —L m , D C A = ^D c m , (4) 

where 

4 

m=l 



\ 



E 

m=l 



D 



\ 
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the SM relevant superpotential becomes to be 

W = fiHA + ^DlD^ + yijLiHiEf + y^QiH^ + yijQiHJI^ 

•A, ; /,/.,/, ; /•;;,: + x'^Q^di + A/jg, /.,/;•, + yfQ % u d Di + y^D^ , 



(6) 



where field decomposition have been generally written as follows, 

L m = c mi Li + c^iTd , D c m = c^DI + c® 4 Dl , (7) 
and the coefficients are 



y%j 2A JT1) ^jCyjjjC n ,4 , Aj mrl C m 4C n j , Ajj'fc X mn kC m iCnj i 

\l - \> r r D \ D — \' r r D n D — V r r D 



(8) 



From the superpotential we see that because of Dirac mass terms of up- type Higgs 
and the four doublet leptons, D C H and the four singlet down-quarks, one of the four lepton 
doublets and one of the down-quarks, namely the fourth doublet lepton H d and the fourth 
singlet down-quark D\ are always heavy, H d is identified as the down-type Higgs. The fourth 
neutrino together with the "neutrino" in H u consists of neutral Higgsinos. After the mass 
terms, the next five terms in Eq. ([6]) are ordinary Yukawa interactions and trilinear lepton 
number (R-parity) violating terms. The other three terms in fl6]) are new which involve the 
D c A ^ H -j field. Two of them also violate lepton numbers. 

Soft SUSY breaking mass terms should be included into the Lagrangian. In addition to 
gaugino masses, they include mass-squared terms of scalars and -B/i-type terms correspond- 
ing to those //-terms in superpotential (J3]), 

-C D M 2 LlL m + MthlK + M%E?Ej + M^Q, + MtfJ?Vf + M*D*D° m 
+M 2 DH D<gD c H + (Bfi m L m h u + B D ^D c m D c H + h.c.) , 

where tildes stand for scalars. We have assumed universality of the mass-squared terms and 
the alignment of the B terms, namely both the mass parameters B and B D do not depend 
on the sub-script m. In terms of three light generations of Eq. ([6]), universality of these soft 
mass terms is easily seen, 

-C D M 2 l\U + M 2 h\h d + MlhlK + M%EfEi + M^Q|q. + M 2 UfU^ + M 2 D D%D° m 
+M 2 DH D%D C H + (Bfxh d h u + B D fx D DlD c H + h.c.) . 

(10) 
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Numerically soft masses M's, -B's and gaugino masses are assumed to be 0(100) GeV. Soft 
trilinear terms corresponding to Eq. ([3]) are 

C D X m mL m L n E^ + X' imn Q t L m D c n + y ij Q l h u U c j + § i:j E°D c H U° + h.c. , (11) 

where the following coupling alignment will be assumed, 

Amni = A mn jm , Kmn = Kmn m , Vij = Vij^O (12) 

with m being the order of soft masses ~ 0(100) GeV. 

Let us look at gauge symmetry breaking. From the Lagrangian, the scalar potential can 
be written down straightforwardly. To get EWSB, one needs a negative determinant of the 
Higgs mass-squared matrix, namely 

(M 2 + /U 2 ) (M 2 + yU 2 ) < \B^\ 2 (13) 

with the ordinary condition M 2 + M\ + 2/i 2 + 2Bfi > 0. This requirement can be realized 
when the renormalization group is considered. M 2 will become negative at the weak scale, 
due to the large top quark Yukawa coupling. Therefore, everything of EWSB here will be 
the same as that in MSSM. The MSSM analysis of EWSB applies here. EWSB in this model 
occurs at the weak scale. Besides Eq. ( TlBl) . correct EWSB also requires 

{Ml + ^ 2 ) (M 2 DH + ^ D2 ) > \B D fi D \ 2 . (14) 

Then the remaining analysis of EWSB is identical to that of MSSM with same Higgs and 
Higgsino spectra. Eq. (fl4l) can be satisfied easily. Carefully thinking of EWSB conditions 
Eqs. (ITBl and (JI4"|) . we see that if \i < fi D , EWSB occurs naturally. 

B. Dark sector 

The dark matter sector Lagrangian is written according to the gauge invariance, 
£dark = (xlePK+giVi+tivi^ + x l 2e - 92 v 2 - giVl - g [v{ X2 + ^aMfa + 0t e -W 02 + x t x j \ gg 

+ (n'XiXs \ee + cX(0i0 2 - fi" 2 ) \ ee + h.c.) . 

(15) 

where \il and //" are mass parameters, g[ and c coupling constants. It is important to note 
that an accidental Z 2 discrete symmetry appears, under which \\ an d X2 fields are odd and 
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all the other fields are even. As we will see, it remains unbroken even after U(l) n breaking 
as well as EWSB. The Z2 symmetry keeps the lightest component of Xi = 1)2) stable. 
That is the dark matter particle in this model. 

When SUSY breaking is considered, soft masses should be included, 

£dark,soft = -Im'^'A 1 ' + m^xtxi + m l 2 X*2X2 + m^fafa + m\ 2 fa 2 fa 
+m 2 x x*x + {B'p!)xiX2 + h.c.) . 

U(l) n gauge symmetry breaks spontaneously when </> lj2 get non-vanishing vacuum expecta- 
tion values (VEVs). From Eqs. (TT5|) and (fl6|) . the relevant scalar potential for and (fa) 
determination is 

U dark = 2g[ 2 ((fa) 2 - (fa) 2 ) 2 + c 2 \(fa)(fa) - // 2 | 2 + ml ((fa) 2 + (fa) 2 ) , (17) 
where it has taken m 2 = m 2 ^ = m 2 2 . Then 

/ m 2 \ 2 

(&> = <&> = <*>= U" 2 -— • (18) 



a 



c 



Vector-like particle masses are all taken to be similar [8j, hence //, //' ~ 1 TeV. It is natural 
to expect that (fa ~ (9(100) GeV. The U(l) n gauge boson j n gets a mass of m ln = 4g[(fa. 
As long as the bosonic fields of SU(2) L doublets Xi,2 do not get VEVs and are heavy enough, 
the Z2 symmetry still remains after U(l) n breaking. 

In addition to Eqs. (IT5|) and ( I16p . the Lagrangian should include a gauge field mixing 
between U(l) n and U(l)y, 

^mixing ^-^n ^Y^tu • (1^) 

This mixing is accompanied with a gaugino mixing because of SUSY, 

£gaugino mixing = 2e(AVd /1 A 1 ' + A^^A 1 ) . (20) 

It is conventional to choose the mixing to be of the size of e ~ 10~ 3 . This mixing makes 
possibly lighter particles in the dark sector, such as fa >2 , decay into MSSM particles. 

What we are interested in is the spectrum of xi an d X2 particles, because they carry 
SM quantum numbers. For the fermions, at the leading order, they form a Dirac particle 

\1> X = ( 1 I with a mass //, 



X2 



£x = W"*x - m'*x*x , (21) 



where = — ig 2 A^r a — ig\A 1 ^ — ig[A'^. Actually \P X is a pseudo-Dirac particle because 
of gauge symmetry breaking. Generally, EWSB splits the neutral and charged components 
of \l/ x , and U(l)„ breaking further splits the two neutral components. In this model, such 
mass splittings are described by the following gauge symmetric dimension 5 and 6 operators, 



C D £ dim - 5 + £ dim - 6 , (22) 



where 



£ dim.5 = ?±( XlHu) ( X2Hd) \ 9Q + ^( XlHd) ( X2Hu) \ ee + h a ? 

£ dim. 6 = ^L MxiHu) ( XlHu) \ ee + ^(^H^Ha) \ ee + h.c. 



(23) 



with A being a cutoff, and ai ~ 0{1) coefficients. A may be considered as the scale of SUSY 
breaking messengers which have been integrated out, the messengers also form complete 
SU(5) representations and do not break the unification of SM gauge couplings. As it will be 
seen, £ dim - 5 splits charged and neutral components, and £ dim - 6 splits the two neutral ones. 

Note that without dimension-5 operators, EWSB itself at the renormalizable level induces 
mass splitting between charged and neutral components. One-loop diagrams with a Z boson 
propagating in inner lines directly give splitting roughly as f^-Mz. It is of the order 0.1 
GeV. For splitting the two neutral parts, we nevertheless need the new higher dimensional 
operators of £ dim - 6 . In such a situation, it is natural to expect that a general £ dim - 5 is also 
there. In our analysis £ dim - 5 is taking as the effective operators which parameterize all the 
EWSB effects. 

ai and a 2 terms simply give rise to mass splitting between the charged and the neutral 
components, 

* , r / x v<2 sin 28 .„ JX 

AM = (a 1 + a 2 ) — P - , (24) 

where v = 246 GeV. With tan/3 between (3 - 10) and A between (10 - 100) TeV, this 
splitting ranges from (0.1 — 1) GeV. Here the positive (ai + a 2 ) case is chosen so that the 
neutral components of Xi, Xi an d X21 are lighter. The charged components Xi an d xt form 
an exact Dirac particle. We can see that mass splitting generated by dimension 5 operators 
is at least as large as that generated by loop diagrams. 

Xi and x 2 are splitted further through the a 3 and a 4 terms. The mass matrix of Xi, X2 
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turns out to be 



M 



( 

1 
\2 



■A } 2 
f 2 sin 2/3 



4A 



a 3 sin 2 (3 
ai + a 2 ) 



' , v 2 sin2/3. 1 \ 

A* ^ — (°i + °2) 

( V ^(<t>) 2 fl 

(^) ^-a 4 cos (5 j 



(25) 



As the off-diagonal elements are much larger than the diagonal ones, the mass eigenvalues are 

, , V 2 Sin 2(3, \ i /■ U \2 (0) /■ -In 2 n\ ii i- 

approximately fi — (ai + 02) ± (-H — (a 3 sin p + a 4 cos p), and the corresponding 

mass eigenstates are approximately 

X / d = (Xi + ^)/v / 2 + 0(^), X^^-^/v^+O^) (26) 

with being the lighter state. Therefore, mass splitting between the two neutral Majorana 
fermions is 

Am= (£) (0>(a 3 sin 2 /3 + a 4 cos 2 (3) . (27) 

This splitting is almost independent of tan/3. Taking (0) ~ 100 GeV, this splitting ranges 
from (1 - 100) MeV. 

The other particles of the dark sector have the following spectrum. As </>i, 02 get VEVs, 
their fermionic partners 0i,2, X which is the fermion of X, and gaugino A 1 ' will get masses. 
It is convenient to change the basis to 



(28) 



$' = (0 X + 2 )/V2 , $ = (0i - 2 )/V2 . 



$' and X form a Dirac particle with a mass of \/2c(0) ~ 100 GeV. The mass matrix of $ 
and A 1 ' is the following by further considering the A 1 ' soft mass, 



M' = 



-lAg'M 
-iAg[ ((f)) m[ 



(29) 



The matrix elements are all (9(100) GeV, so the mass eigenstates A" and N' are of the same 
mass scale, with A" the lighter one. 

The scalars xi,2 and 0i )2 are heavy ~ 100 GeV —1 TeV due to soft SUSY breaking. Notice 
that the singlet x has no gauge couplings and can have a vanishing soft mass in the case of 
gauge mediated SUSY breaking. In this case the boson x will have a mass ^2c(0) ~ 100 
GeV which is degenerate to the corresponding fermion. 
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It is convenient to express the fermions in the 4-component form: 

*-=| ,%=[ X ' d ) ^-=( Xl - I ■ (30) 



Xd ) \ X d ) \ X2^ 

In terms of all above mass eigenstates, the dark sector Lagrangian relevant to dark matter 
annihilation can be expressed as 

/COS(7\y 

-g 2 sm6 w A^^-) (31) 

It is seen that the dark matter particle Xd mainly scatter inelastically via gauge interactions. 
Note that there are still small diagonal x^-Xd-gauge boson couplings which are about e' ~ 
Am/fj,' as can be seen from Eq. (1261) . 

C. UV-completion 

Up to now, our model is a TeV effective theory which does not include particles much 
heavier than TeV. The cut off A may be, as we have mentioned, understood as a result of 
integrating out SUSY breaking messengers in the gauge mediated SUSY breaking scenario. 
Here we present an UV completion model which reproduces our effective theory. 

The SUSY breaking messengers have the following quantum numbers of 
SU(2) L xU(l) y xSU(3) c xU(l) n , 

t?(2, -1,1,-1), r/(l, §, 3, -1) , 0, 1, -1) , 

(32) 

fj{2, 1, 1, 1) , fj'{l, -I 3, 1).«(1, 0,1,1), 
The SUSY breaking spurion field S couples with messengers 

y^spurion = S [Tjfj + 7/7/ + kk) , (33) 

where S get a VEV, 

(s) = A + eeF. (34) 
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The soft masses are the same as those of normal gauge mediation. For dark sector, 

a[ F 



< 2 - & (35) 
The U(l)„ messengers k and k have Yukawa couplings with Higgs and dark sector, 

Wmessenger = hXlH u R + b 2 X2H d K + b 3 (j) 2 KK + b^KK . (36) 

These couplings give rise to the tij couplings on tree level by integrating out k and k, 

cti = bib 2 

a 3 = hb\ (37) 

04 = 6462 • 

Here the lack of a 2 will not change mass splitting result of the dark sector significantly 
because a\ is non-zero. There could be models in which a 2 does not vanish. 

The Yukawa couplings bj cause small mixing between messengers and dark sector, which 
makes messengers decay into dark sector. After EW and U(l) n breaking, a Z 2 symmetry 
remains in messenger and dark sector, which is the same X 2 that appears only in dark sector 
when messengers are integrated out. 

III. ANALYSIS 

A. Dark matter relic density 

From discussions of the last section and Eq. ( 13TI) . we see that the only stable particle in 
this model is Xd, which is dark matter. 

In relic density calculation, coannihilation of all the four components of Xi an d y ? should 
be considered. We use the program micrOMEGAs to calculate relic density [20I [2 ll] {22 1. As 
long as Xd has the correct relic density, the dark matter mass fi' can be determined by the 
coupling g[. 

From FIG [1] we can see that the dark matter mass is between jj! = (1.1 ~ 1.5) TeV when 
g[ ranges between to g 2 . This means that taking g[ ~ g 2 will not deviate too much from 
normal WIMP expectation. 

11 



DM Mass 



1500 



1450 - 
1400 - 



< 1300 - 

3. 

1250 - 

1200 - 

1150 - 




0.2 0.4 0.6 0.8 1 

gi'/g 2 



FIG. 1: The relation between dark matter mass // and g[/g2- 
B. Detection of Dark matter 



Direct detection experiments for dark matter have given strict upper limits on cross 
sections of its scattering with nucleus. As Xd mainly interacts with Z boson inelastically, 
it is possible to suppress the scattering cross section with nucleus. It has been known that 
if the mass splitting Am between Xd and x'd * s zero, the x<r n ucleon spin-independent cross 



section would be about 10 38 cm 2 



llj while XenonlOO gives the upper limit of the cross 



section to be 10 -44 cm 2 for ~ 1 TeV dark matter [12j. In our model, Am ranges from (1-100) 
MeV depending on the cutoff A and is larger than the possible kinematic energy of dark 
matter of a maximum speed of 600 km/s. Therefore dark matter will not scatter through 
vector interaction with Xe nuclei on tree level. 

The axial-vector coupling of Xd with Z still needs consideration. This coupling is elastical 
and of the size O(^jr), which is between 10~ 4 to 10~ 6 . It will bring to spin-dependent cross 
section with Xe. The spin-dependent tree level cross section for WIMP elastically scattering 



at zero momentum transfer is 



13|: 



°sd = 32(J ^ 1] G 2 F M?\a p (S p ) + a n (S n )\ 2 (38) 

7Tj 

where J is nuclear spin, M r is the reduced mass of dark matter and target nuclear, a P)Tl 
are effective proton and neutron couplings, (S p>n ) are spin expectation values of proton and 
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neutron in the nuclear. The a PiTl for \d is of the mag nitude The Xd-Xe spin-dependent 
cross section is (10~ 45 — 10~ 41 ) cm 2 for Xe 129 and (10~ 46 — 1CT 42 ) cm 2 for Xe 131 . These cross 
sections are far beyond the detection capability of XenonlOO or Xenon It. 

At one-loop level, the dark matter-nuclear cross section is essentially the same as that 



discussed in Ref. 



14(: the one-loop x^-nucleon spin-independent cross section is about 10 



-48 



cm 2 , which is also too small to be detected. 

Dark matter in our model cannot produce positron excess observed by cosmic-ray experi 



ments such as PAMELA and Fermi-LAT 
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The reason is that there is no light particle 



to provide enough Sommerfeld enhancement for dark matter annihilation in our model. As 
a result, we have to consider the observed positron excess as an astrophysical phenomenon. 



C. Collider phenomenology 

Experimental constraints should be considered. In this model, particles beyond those 
of MSSM, which also carry color or electric charges are D c Al D c H and XiiXt- Each pair 
forms an exact Dirac particle. The former is down-type quark like, and the latter charged 
lepton like. The direct experimental search at LEP requires that they should be heavier 
than 100 GeV, that at Tevatron for down-type heavy quarks requires they are heavier than 
372 GeV 15j, and at LHC the current limit is 675 GeV for the down-type heavy quark 
16l | . These results can be simply satisfied if fi D is larger than 675 GeV and y! larger than 
100 GeV. The electroweak precision measurements generally have weak constraints on this 
model, because these extra matters are vector-like which have contribution in the form of 
l/fi D ^ 2 as expected from the decoupling theorem. The effect of the extra matters can be 
small enough < (m t /fi D ^) 2 ~ (1 - 10)% if we take fi D ® ~ 500 GeV -1 TeV. Noting that 
these direct search limits are obtained with assumption of single decay channel dominant, 
we will take that fi D = 500 GeV and // ~ 1.1 TeV in numerical illustration. 

There are constraints comin g fr om the unitarity of the 3x3 CKM quark mixing ma- 



trix of three chiral generations 15J. This unitarity is consistent with current data within 
experimental errors. In this model, extra down-type quarks mix with ordinary three chiral 
down-type quarks, which necessarily break the unitarity of the CKM mixing matrix. Uni- 
tarity violation is about {ran/ y D ) 2 where the mass parameter = —yfv cos/3/a/2 as seen 
from Eq. ([6]). This jj, D dependence is generally expected in the case of extra vector-like 
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quarks. Hierarchical or small mixing masses can easily make the CKM matrix approx- 
imately unitary within errors. For an example, (mu/fi ) 2 < 10~ 3 . Assuming only the 
third generation mixes with extra quarks, the constraint is still loose, (m 34 //i D ) 2 < 0.39. 
The quantity m 34 is at most about m t . This gives that the parameter fi D > 280 GeV. It 
is easy to see that there are new phases in fermion mixing matrices. However, these new 
matrix elements are of order of (m t //i D ) 2 at most. So new CP violation effects are generally 
suppressed. 

Decay signals of these new particles can be easily identified. From trilinear Yukawa 
interactions given in Eq. (jH]), it is seen that D c ^D c n decay into SM first three generation 
matters. Denoting the new Dirac quark in (D^, D^) as ^p, decays of have following 
results, 



16n 



D\2 



(39) 



Taking relevant Yukawa coefficients yf's ~ 10 _1 — 10~ 2 , the decay rate in Eq. (1391) is 
T ~ 5 — 500 MeV. Taking EWSB into consideration, d mixes with SM fermions. We see 
that the decay — > i W + occurs via the SU(2)i gauge interaction at the level of 0{m/ fi D ), 



G F m 2 w \n D \\V, 




where V35 is the mixing element, and the phase space factors were given in Refs. [171 ] . Taking 
m/n D ~ 1/3, the T is about 1 GeV. In this decay process, the top quark further decays 
into three quarks with one of them a bottom, and the W can decay into a single charged 
lepton with a neutrino. Taking this process as the main decay channel, in the case of ^ d 
pair production, the signal can be searched in the events of 2 charged leptons (electron or 
muon) and 6 jets with two of them b-jets, and large missing energies. 

Denoting the heavy charged Dirac lepton (x\,xt) as ^ s m ass is only AM ~ (0.1 - 
1) GeV above the dark matter. It can decay into or ^>' d together with a pion or with a 
charged lepton (electron or muon) and a neutrino. In the limit of AM << //, the typical 
decay width of is 

r(¥_ e" * e ) ~ . (41) 
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This lifetime is about 10~ 12 — 10~ 7 s. For a larger AM, decay rapidly into a charged 
lepton and missing energy, in which situation detection is difficult. For AM < 0.3 GeV, 
is long-lived and can leave tracks in detectors. As is always produced in pairs, the signal 
should be easy to identify. 

The new quark can be produced at LHC via the gluon fusion process, g a — > ^d^d- 
The production mechanism is essentially the same as that of the top quark 18| with an 
estimated cross section ~ hundreds fb by taking fi D ~ 500 GeV and y/s = 14 TeV. For the 
new charged lepton, the Drell-Yan process is the main production mechanism. The cross 
section is estimated to be a few fb which means a few events in one year at most • 

IV. DISCUSSION 

Finally we will discuss some physical aspects related to this model. First of all, it is 
important to discuss GUT. We have kept gauge coupling unification at the high scale, but a 
real GUT model with a simple gauge group has not been given, it is still far from our reach. 
From the particle content and assumptions of this model, it is doubtful if such a GUT model 
exists. We wonder if a simple gauge group is really necessary for the so-called unification. In 
fact, in certain string models, unification is achieved without requiring a GUT gauge group, 



but the gauge coupling constants are unified at the scale of 10 16 GeV or so [231]. It would 
be nice if this model can be reconstructed as a string model. 

The assumption of baryon number conservation may have a better looking. It has been 
shown that it can be replaced by that of a Z 3 discrete symmetry which is called baryon 



parity : 



241 ]. As we know that any global symmetry is not favored from the point of view 
of the quantum gravity, because black holes violate such symmetries. However, the baryon 
parity can be considered as a result of gauge symmetry breaking [25( . 

The two U(l) gauge interactions have a mixing. With such a mixing which is symmetry 
allowed, it is guaranteed that dark matter is composed of Xd only. The other particles of 
the dark sector can decay via the mixing even if they are lighter than Xd- This is a simple 
scenario of dark matter in this model, although a complicated one with vanishingly small 
e is possible. Without the mixing, the lightest particle of the dark sector other than xi,2 
would be stable and contribute to the relic density of dark matter. In that case, dark matter 
would be composed of multiple components. 
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Our dark matter model has nothing to do with the so-called indirect indication of the 



dark matter by astrophysical observation of ATIC [26] , PAMELA [27] and FermiLAT 



28|. It 



looks that our model has a potential to accommodate their observation. They have observed 



an access of cosmic electrons and positrons. Although Fermi-LAT experiment [28j does not 
support ATIC, it agrees with PAMELA. The observation has inspired a lot of theoretical 
re-consideration about the WIMP dark matter jf], 30] . Arkani-Hamed et al. [s| have 
proposed a scenario to understand all experiments on the dark matter in a natural way. 
WIMPs should have new interaction mediated by a light, GeV scale particle which enhances 
their current annihilation cross section via the Sommerfeld mechanism. This GeV particle is 
supposed to be the main annihilation product, and because of its lightness, it finally decays 
into leptons only. There are many theoretical models realizing this new scenario |29j . There 
is a point which is similar to our case. In our model, the WIMP also has new interaction 
which, however, has a typical energy scale of hundreds GeV. If we had made some tuning 
on the mass parameter //" to reduce the new interaction scale to be GeV, this model would 
realize most part of the scenario of Arkani-Hamed et al. But because of the dark matter 
structure of this model, the Sommerfeld enhancement factor would be about 10 4 which has 
been ruled out [311 ] . 



V. SUMMARY 



Within the framework of R-parity violation, we have studied the dark matter problem 
with the constraint of gauge coupling constant unification. The WIMP dark matter is 
contained in a new vector-like SU(2)l doublet which possesses a new U(l) gauge symmetry. 
The Higgs particles are included in a 5 © 5 representation of SU(5). Instead of R-parity, 
baryon number conservation is assumed. In this model, the dark matter particle is stable 
as a result of the gauge invariance. An accidental discrete Z<i symmetry remains after 
spontaneous gauge symmetry breaking, which makes the dark matter particle stable. 

Main results of this model is the following. The mass of the dark matter particle is 
(1.1 — 1.5) TeV. The dark matter and nucleus interaction has small cross sections which 
are consistent with the current dark matter direct detection experiments like XenonlOO. In 
addition to the particle content of MSSM, we have a new down-type Dirac quark and a 
new Dirac charged lepton with masses about (500 — 1000) GeV. These charged particles, 
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especially the new quark can be produced at the LHC and hopefully be observed. 
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